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(N ! ABSTRACT 

Context. The relationship between the emission-line equivalent width and the continuum luminosity, called the Baldwin effect, 



> 
ON 



is important for studies of the physics of the broad-line region of AGNs. Some recent studies have revealed the non-constant 
slope of the intrinsic Baldwin effect for several Seyfert 1 galaxies. 

Aims. Using the archived ultraviolet spectra obtained by IUE, HST, and HUT in 1978-2002, we investigated the intrinsic 
, Civ Baldwin effect of the well-studied Seyfert 1 galaxy NGC 4151. Both its continuum flux and Civ emission-line flux varied 
about two orders of magnitude in more than two decades, making it one of the best targets for studying the slope variations 
\£i caused by the Baldwin effect. 

f*"*) Methods. We fitted the Civ line profile of the 490 archived UV spectra of NGC 4151 with a similar model consisting of a few 
Gaussian components and derived the slope in the log-log plot for the total flux of C IV emission line against the UV continuum 
flux in different observation epochs. 
I > Results. We found that the slope is not constant for NGC 4151, as it varies from 0.58 in the highest flux epoch to 0.83 in the 
O lowest flux epoch. The slope evidently decreases as the continuum flux increases, which reinforces the previous findings of the 
+3 ■ non-constant slope in the H/3 Baldwin effect of NGC 5548 and the C IV Baldwin effect of Fairall 9. 

ffi ' Conclusions. Our result suggests that such a non-constant slope may not be unusual for AGNs. Its physical origin is probably 
. . related to the different non-linear responses of the emission line to the variable ionizing continuum caused by the different 
■ accretion modes at different luminosity levels. We briefly discuss the effects of various absorption components in the C IV line 
profile of NGC 4151 and argue that the slope variation is not driven mainly by the absorption effect. Intensive, long-term, and 
high-resolution spectral observations of strongly variable AGNs are needed in the future to confirm our result. 
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1. Introduction 

The relation between the continuum and emission-line lu- 
minosities is important for our understanding of AGN 
physics. Such a relationship, known as the Baldwin effect, 
states that the broad emission-line equivalent width (EW) 
decreases as the continuum luminosity increases (Baldwin 
1977; for a review see Osmer & Shields 1999). The Baldwin 
effect was first discovered for the C iv A1549 line (Baldwin 
1977) and was confirmed later for most of the other strong 
UV emission lines (Kinney, Rivolo & Koratkar 1990), as 
well for the optical hydrogen Balmer lines (Gilbert & 
Peterson 2003; Goad, Korista & Knigge 2004). 

The physical origin of the Baldwin effect is still un- 
clear, but is probably related to the luminosity-dependent 
continuum variations. According to a theoretical work by 
Korista, Baldwin & Ferland (1998), the ionizing contin- 
uum shape is softened if the luminosity increases. The 
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reduction in the ionizing photons at a given UV/optical 
luminosity leads to smaller equivalent widths of the broad 
emission lines. Indeed, some accretion disk models can also 
explain the continuum softening naturally in the case of 
increasing luminosity (Netzer, Laor & Gondhalekar 1992). 
However, other factors, such as the selection effect, vari- 
ability, and light travel-time effects, may also account 
for some parts of or the scatters of the Baldwin effects 
(Jones & Jones 1980; Murdoch 1983; Krolik et al. 1991; 
Pogge & Peterson 1992; Peterson et al. 2002). Recent stud- 
ies, including Baskin & Laor (2004) and Bachev et al. 
(2004), have revealed the possible correlations between the 
C iv EW and some parameters defining the Eigenvector 1 
of AGNs (Boroson & Green 1992). Shang et al. (2003) also 
show that the Eigenvector 1 can contribute to the scatters 
of the Baldwin effect. 

The Baldwin effect can be expressed by the simple for- 
mula EW oc L@, where L c is the continuum luminosity, 
or alternatively by Lu nc oc L™, where Lu ne is the emis- 
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sion line luminosity (evidently a = 1 + /?). There are 
two kinds of Baldwin effects treated in the literature. One 
is the global Baldwin effect, which is obtained from the 
single-epoch observations for an ensemble of AGNs. The 
other is the intrinsic Baldwin effect, which represents the 
line-continuum relation in a single variable AGN. For the 
global Baldwin effect, the (5 value was found to be —0.17 
and —0.12 for Civ and Ly a lines, respectively (Kinney 
et al. 1990; Pogge & Peterson 1992). Some studies also 
indicate steeper slopes of the global Baldwin effect for the 
relatively high ionization lines than those found for the low 
ionization ones (Wu et al. 1983; Kinney et al. 1987, 1990; 
Baldwin et al. 1989; Zheng et al. 1997). The flattening 
of the slope of the global Baldwin effect at lower contin- 
uum luminosities was mentioned previously by Osmer & 
Shields (1999) as a second-order effect, and was also no- 
ticed recently by Baskin & Laor (2004) from a study of 81 
BQS quasars. 

The intrinsic Baldwin effect (usually expressed as 
fline K where fu ne and f c are the emission line 
and continuum fluxes, respectively), has been studied only 
for several strongly variable AGNs. The derived a value 
varies from 0.1 to 0.6 for different AGNs, with a mean 
value around 0.4 (Kinney et al. 1990; Pogge & Peterson 
1992). Recently, Gilbert & Peterson (2003) and Goad et 
al. (2004) studied the intrinsic Baldwin effect for the best- 
studied Seyfert 1 galaxy NGC 5548 using the data of 13- 
year-observations and found that the a value for the H/3 
line varies in a range from 0.4 to 1. Goad et al. (2004) 
find that the a value decreases as the continuum flux in- 
creases, suggesting the slope of the intrinsic Baldwin effect 
is not constant. Wamstekcr & Colina (1986) and Osmer & 
Shields (1999) also notice the slope change of the intrinsic 
C iv Baldwin effect for another Seyfert 1 galaxy, Fairall 9, 
which clearly shows the flattening of the EW value as the 
continuum flux decreases. Interestingly, the trend of the 
slope change in the intrinsic Baldwin effect seems consis- 
tent with that of the global Baldwin effect. However, this 
may not indicate that the physics behind them is totally 
the same. For example, either the different metallicity or 
black-hole mass has been invoked as possible origin of at 
least part of the global Baldwin effect (Warner, Hamann 
& Dietrich 2004), but this is clearly not relevant to the 
intrinsic Baldwin effect of a single AGN where these pa- 
rameters are fixed. 

To further test the non-constant slope of the intrin- 
sic Baldwin effect, in this paper we investigate the rela- 
tionship between the C IV emission line flux and the UV 
continuum flux for another well-known nearby Seyfert 1 
galaxy NGC 4151 that has been extensively observed in 
the UV band by the International Ultraviolet Explorer 
(IUE), Hubble Space Telescope (HST), and Hopkins 
Ultraviolet Telescope (HUT) in the past three decades 
(Boksenberg et al. 1978; Clavel et al. 1987, 1990; Ulrich 
et al. 1991; Crenshaw et al. 1996, 2000; Kriss et al. 1992, 
1995; Weymann et al. 1997; Kraemer et al. 2001; see Ulrich 
2000 for a detailed review) . Both the UV emission- line flux 
and continuum flux of NGC 4151 varied about two orders 



of magnitude in this long observation period, making it 
one of the best targets for studying the variation in the 
slope of the Baldwin effect. 

In Sect. 2 we present the data analysis of the UV spec- 
tra from IUE, HST, & HUT. In Sect. 3 we show our result 
for the varying slope of the intrinsic Civ Baldwin effect. 
Finally we give our conclusions and briefly discuss our re- 
sult in Sect. 4. 



2. Data analysis of the archived UV spectra 

To investigate the slope variation of the intrinsic 
Civ Baldwin effect in NGC 4151, we measured the flux 
of the C iv emission-line and the flux of UV continuum 
observed in a long period. 

2.1. The data set 

As a nearby (cz = 995 km s _1 ) bright Seyfert 1 
galaxy, NGC 4151 has been extensively observed by IUE 
(Boksenberg et al. 1978; Clavel et al. 1987, 1990; Ulrich 
1996; Ulrich et al. 1991; Crenshaw et al. 1996; Edclson 
et al. 1996). However, different results have been found 
for the correlation between the UV continuum and the 
emission-line flux variations. For example, from the obser- 
vations in the 1993 IUE campaign Crenshaw et al. (1996) 
found that the emission-line light curve does not correlate 
well with the continuum variation over the short duration 
of observations. But from the long-term light curves over 
two decades, the response of the emission-line flux with 
the continuum variation was clearly observed (Ulrich et 
al. 1991; Ulrich 2000). In this paper, we re-investigate this 
correlation using the 468 archived low-resolution large- 
aperture spectra of IUE/SWP from 1978 to 1996 (JD 
2443722.59 to JD 2450243.40) where the Civ emission 
line was clearly detected. 

Furthermore, NGC 4151 was observed by HUT at a 
resolution of 3 A in December 1990 and March 1995 (Kriss 
et al. 1992, 1995), and the 8 archived HUT spectra are 
available. These observation were done in a time period 
that partly overlapped with that of the IUE observations. 
Therefore, the HUT spectra can be used to test the mea- 
surement accuracy of the IUE data. 

NGC 4151 has also been frequently observed by HST. 
The 15 HST/STIS archived spectra in 1998-2002 are 
available and used in our study. However, the archived 
HST/GHRS spectra were skipped because they do not 
cover the whole red part of the C IV line profile (Weymann 
et al. 1997). In all of the 15 HST/STIS spectra, the whole 
C iv line profile is shown clearly. 

In total, we collected 490 archived UV spectra of NGC 
4151 with a wavelength range from about 1300 A to 1800 A 
observed by IUE, HUT, & HST in the period of 1978-2002. 
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2.2. The continuum and Civ emission-line flux 
determinations 

Low-resolution spectra of NGC 4151 taken with IUE and 
HUT clearly show broad emission and remarkable cen- 
tral absorption in the C iv line profile. However, high- 
resolution spectra from HST have revealed many more 
other subtle emission and absorption features (Wcymann 
et al. 1997; Crenshaw et al. 2000; Kracmcr ct al. 2001). 
Therefore, the Civ line profile of NGC 4151 is very com- 
plex and needs to be treated with caution. Because 468 
of the 490 archived spectra in our study were obtained by 
IUE, the lower resolution of IUE spectra does not allow a 
more accurate fitting of the subtle emission and absorp- 
tion features. In order to avoid the systematic uncertain- 
ties caused by different fitting models, in this study we 
simply adopt the same spectral-fitting model for the HST 
and HUT spectra as for the IUE spectra. 

The accurate measurement of the continuum flux is 
important. Usually, the local continuum fitting with a 
power-law or a straight line can produce reasonable mea- 
surement if the fitting is made with caution. Here we use a 
straight line to fit the local continuum in the selected con- 
tinuum windows, namely 1260-1290 A, 1420-1460 A, and 
1805-1835 A. The results are satisfactory. In fitting the 
continuum, we find that the C n A1334A line has rela- 
tively strong effects on the continuum flux measurement at 
1350 A. Therefore, we adopt the continuum flux at 1440 A 
for our study of the line-continuum relation of NGC 4151 
because the lines around 1440 A, like O iv] A1402A and 
N iv] A1486A, are relatively weak. The continuum flux 
at 1440 A was taken as the weighted mean over the wave- 
length range 1420- 1460 A, and the flux uncertainty is esti- 
mated from its standard deviation. No correction is made 
for the extinction, as it is negligible for NGC 4151 (Kriss 
et al. 1992). 

In the IUE spectra of NGC 4151, the Civ AA1548.2, 
1550.8 lines arc blended with a few variable components 
covering a wavelength range from 1480 A to 1700 A ( sec 
Clavel et al. 1987 for the details). To reduce the contam- 
ination from other lines, we fit the continuum-subtracted 
spectra within the wavelength range from 1450 A to 
1720 A simultaneously with ten individual Gaussian com- 
ponents. Five lines like N iv] A1486A, LI A1518A, L'2 
A1576 A, L2 A1594A, O ill] A1663A are fitted with a sin- 
gle Gaussian component. The C iv line is fitted with three 
Gaussian components: a broad emission component, a nar- 
row emission component, and a narrow absorption com- 
ponent. The H II A1640A emission line is fitted with two 
Gaussian components including a narrow one and a broad 
one. Two narrow and variable satellite emission lines (LI 
and L2) were first noticed by Ulrich ct al. (1985) when 
NGC 4151 was at a low continuum state. Using the high- 
resolution HST/STIS spectra of NGC 4151, Crenshaw et 
al. (2000) indicate that these satellite emission lines (LI, 
L'2, and L2) shown in the IUE spectra are not in fact 
separate lines, but rather the residual emission that ap- 
pears isolated due to the low-ionization absorption lines 



(primarily Si n A1526.7A, Si II* A1533.4A, and Fe n mul- 
tiples) overlaying the large part of wings of C iv. Assuming 
them to be the satellite emission lines can slightly under- 
estimate the total emission line flux of C iv in the low- 
continuum case. The effects of such a simplification on 
our result will be addressed in Sect. 3. 

The nonlinear least-square Levenberg-Marquart min- 
imization method is used in our spectral fitting. Unlike 
Clavel et al. (1987) who fixed the line widths of some 
Gaussian components, we allow the fitting parameters for 
all components, such as the line central wavelength, veloc- 
ity dispersion, and the area under a line, to vary at certain 
ranges. The uncertainties of these parameters are obtained 
using a method similar to the one adopted in Clavel et 
al. (1991). Three examples of our multi-component fitting 
of the spectra taken at the low, moderate, and high flux 
states with IUE, HST, & HUT, respectively are shown in 
Fig. 1. The fittings are satisfactory in general. However, 
sometimes the red peak around 1540-1560 A is too sharp 
to be fitted perfectly, especially for the spectra taken by 
HST. This leads to a slight underestimation of the total 
C iv emission line flux but will not affect our result signif- 
icantly. 

In Fig. 2, we show the light curves of the continuum 
flux at 1440 A, the total Civ line flux (the sum of the 
modeled broad emission, narrow emission, and absorption 
components), and the broad emission- line component flux 
from the IUE, HUT, & HST observations in 1978-2002. A 
huge broad peak in 1991-1998 and several sub-peaks are 
clearly shown in all three light curves. The overall response 
of the C iv emission line with the continuum variation is 
clear in these long-term light curves. This also confirms 
that the emission-line response with the continuum cannot 
be revealed by the short duration observations (Crenshaw 
et al. 1996) but can be clearly detected by the long-term 
spectral monitorings. 

3. The intrinsic Baldwin effect of NGC 4151 

Using the measured flux data for both C iv emission line 
and UV continuum, we can study the slope variation in 
the C iv Baldwin effect. 

3.1. The non-constant slope 

The Baldwin effect represents the dependence of the 
equivalent width (EW) of Civ emission line on the UV 
continuum flux. Figure 3 shows the total C iv emission- 
line flux and the line EW against the UV continuum flux 
at 1440 A for NGC 4151. The data from HUT & HST are 
fully consistent with those from IUE, clearly suggesting 
a non-constant slope in the log-log plots. We notice that 
the trend of the curvature in the upper panel of Fig. 3 is 
similar to that found for the intrinsic H/3 Baldwin effect 
of NGC 5548 (see Fig. 4 of Goad et al. 2004) and Fairall 
9 (see Fig. 2 of Wamstcker & Colina 1986), which also 
show the steepening of the slope at the lower flux state 
and the flattening of the slope at the higher flux state. 
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Fig. 1. Sample spectra of NGC 4151 from IUE (lower 
panel), HST (middle panel), & HUT (upper panel) ob- 
servations and fitting results around the 1450-1720A re- 
gion. The three spectra (from bottom to top) are for the 
low, moderate, and high flux states, respectively. The dot- 
ted line represents the observed spectrum after the contin- 
uum subtraction. The thin solid lines represent the best-fit 
spectral components (see text for details). The thick solid 
line denotes the sum of these components. The residual of 
the fit is also shown in the lower part of each panel. Flux 
is in units of 10 -13 erg s -1 cm _2 A _1 . 



In the EW-continuum plot (the lower panel of Fig. 3), 
the flattening of the slope in the low continuum state is 
also consistent with the trend noticed by Osmer & Shields 
(1999) for Fairall 9. This suggests that the non-constant 
slope of the Baldwin relationship may be usual for AGNs. 

To demonstrate the non-constant slope of the Baldwin 
effect for NGC 4151 more clearly, we divide the obser- 
vations into 4 epochs. Epoch 1 (JD2443722-JD2445654) 
covers the IUE observations in 1978-1983, which shows 
at least two sub-peaks in the light curve. Epoch 2 
(JD2445766-JD2447631) covers the IUE observations in 
1984-1989, which represents the lowest flux state in all 
the UV observations. Epoch 3 (JD2447948-JD2450243) 
covers the IUE observations in 1990-1996 and the HUT 
observations in 1990 and 1995. The number of spectra in 
this epoch is significantly larger than others since it in- 
cludes 203 spectra obtained in a one-month intensive IUE 
campaign in 1993, which covers the highest flux state of 
NGC 4151. Epoch 4 covers the HST observations in 1998- 
2002. Although only 15 HST/STIS spectra are available, 



.3 4 
2 




a 400 

O 
-t-j 

§ 200 



100 




l_ I 

f 


X 



o o: 

1 n I* 1 T 1 


-8 J 1 

Jl g S x a 

7% i x | i x i x 


X 
X 

x o - 

@ Q e" 

1 |°l ? u 


: 1 I 1 1 I □ 

- 1 s 

# ill i 1 1 


x : 

X 

x : 
o £L 



O : 
i 1 Oi®m i- 



2000 4000 6000 8000 
Julian date(-2443000) 



Fig. 2. Light curves of the continuum flux at 1440 A, the 
total C IV flux and the broad emission line component 
of NGC 4151 in 1978-2002. The crosses, squares, and 
circles represent the IUE, HUT, and HST data, respec- 
tively. The average uncertainty was shown in the upper- 
left part of each panel. The continuum flux is in units of 
10~ 13 erg s _1 cm -2 A and the Civ line is in units of 

io- 13 



erg s 1 cm 2 . 



we can still clearly see that both the continuum flux and 
C iv emission line flux fade away in this epoch from the 
major peak. The numbers of spectra, the average UV con- 
tinuum flux, and the estimated slope of the Baldwin effect 
in these 4 epochs are listed in Table 1 . The variation in the 
slope a with the UV continuum flux is shown clearly in 
Fig. 4. The slope a varies from 0.58 in epoch 3 (the high- 
est flux case) to 0.83 in epoch 2 (the lowest flux case). The 
result from the HST observations (a = 0.78 in epoch 4) 
confirms the trend obtained from the IUE and HUT ob- 
servations in the other 3 epochs. A smaller slope (3 value 
(—0.17) in epoch 3 is consistent with the almost constant 
EW value at the lower flux state noticed for Fairall 9 by 
Wamsteker & Colina (1986) and Osmer & Shields (1999). 
The slope a values for NGC 4151 are also well within the 
range for the broad H^ line of NGC 5548 found by Goad 
et al. (2004), who indicated that the slope a varies from 
0.4 to 1 in the 13-year observations of NGC 5548. 

3.2. The effects of absorption features 

The non-constant slope of the intrinsic Baldwin effect has 
been found so far for NGC 4151 (this paper), Fairall 9 
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Table 1. The average continuum flux and slope of the Baldwin effect in four different observation epochs. The standard 
deviations of these two values are also given. The flux is in units of 10~ 13 erg s _1 cm _2 A _1 . 



Epoch JD Time No. of Spectra F ° a(F ■ ) a 8 a 

1440A v 1440V 

1 JD2443722 - JD2445654 75(IUE) 1.077 0.574 0.721 -0.279 0.025 

2 JD2445766 - JD2447631 85(IUE) 0.490 0.307 0.833 -0.167 0.059 

3 JD2447948 - JD2450243 308(IUE) + 7(HUT) 3.607 1.294 0.581 -0.419 0.025 

4 JD2450855 - JD2452403 15(HST) 0.844 0.754 0.777 -0.223 0.056 




-14 -13 -12 

L Og F 1440A 

Fig. 3. Correlations of the total C iv line flux and the 
equivalent width (EW) with the continuum mix at 1440 A 
for NGC 4151. Symbols have the same meanings as in 
Fig. 2. The average uncertainty is shown in the upper-left 
part of each panel. The continuum flux, line flux, and EW 
are in units of erg s _1 cm~ 2 A , erg s _1 cm~ 2 , and A, 
respectively. The curvature of the Baldwin relationship is 
clearly shown in both panels. The results from HUT & 
HST arc fully consistent with those from IUE. 



(Wamstekcr & Colina 1986; Osmer & Shields 1999), and 
NGC 5548 (Goad et al. 2004). However, unlike Fairall 9 
and NGC 5548, which show no absorptions and weak ones 
in their UV spectra, respectively (Crenshaw et al. 1999), 
NGC 4151 shows much more remarkable absorption fea- 
tures, especially in the C IV line profile (Boksenbcrg et al. 
1978; see also Fig. 1). These absorption components also 
vary with the UV continuum. From Fig. 1 we see that the 
central absorption is stronger when the continuum flux 
is higher. This absorption feature has been found to ac- 
tually combine various absorption lines with different ion- 
ization potentials (Kriss et al. 1992; Weymann et al. 1997; 
Crenshaw et al. 2000; Kraemer et al. 2005). Considering 
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Fig. 4. Variation of the slope a with the continuum flux 
for NGC 4151 in 4 observation epochs. The epoch num- 
ber is shown for each point. The epoch 4 is for the HST 
observations in 1998 2002. Sec Table 1 for details. 



the absorption feature, a natural question then arises: is 
the non-constant slope of the C iv Baldwin effect in NGC 
4151 due to the variation in absorption at different con- 
tinuum flux levels? To answer this question, we investi- 
gated the variations in the broad emission, narrow emis- 
sion, and absorption components of C IV line. Figure 5 
shows the variations in these components with the UV 
continuum flux. We can see that the flux values of all 
these three components decrease when the continuum flux 
decreases. The variation in the absorption component is 
very similar to that of the narrow emission-line compo- 
nent. Moreover, we noticed that in the higher continuum 
flux state, the broad-line component always dominates the 
total C IV flux. In the lower continuum flux state, however, 
the narrow emission- line component becomes dominant. 
Both the broad-line component and the central absorp- 
tion components become relatively weak. If we consider 
the Baldwin effect of the broad emission-line component 
alone, we also can obtain a non-constant slope with a vary- 
ing from 0.74±0.02 in the highest flux epoch to 1.33±0.12 
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Fig. 5. Variations in the broad, narrow, and absorption 
components of Civline with the UV continuum flux at 
1440 A for NGC 4151. Symbols have the same meanings 
as in Fig. 2. The average errors are shown in the upper- left 
part of each panel. The continuum flux and line flux are in 
units of erg s _1 cm _2 A _1 and erg s" 1 cm^ 2 respectively. 



in the lowest flux epoch, which clearly indicates that the 
non-constant slope may not be dominated by the absorp- 
tions. 

As mentioned in Sect. 2, besides the remarkable cen- 
tral absorption component, there are also many subtle 
absorption features in the wings of C iv profile of NGC 
4151, as revealed by the high-resolution HST spectra 
(Weymann et al. 1997; Crenshaw ct al. 2000; Kraemer et 
al. 2001). These complex features make it more difficult 
to extract the information of the unabsorbed Civ flux. 
In our spectral fitting, we only use a Gaussian compo- 
nent to represent the central absorption of C iv and treat 
the absorption-induced 'emission' features in the wings 
of Civ as satellite emission lines (LI, L'2, and L2). As 
these features are in fact caused by the low-ionization ab- 
sorption lines overlaying the large part of wings of C iv, 
our simplification would underestimate the total flux of 
Civ. As indicated by Crenshaw et al. (2000), these low- 
ionization absorption features are weak in the high-flux 
state and become prominent when the UV continuum 
drops to a low state. Therefore, the simplification of our 
spectral fitting is fine when NGC 4151 is bright (e.g. the 
continuum flux larger than 10~ 13 crg s _1 cm -2 A -1 , as 
in epochs 1 and 3). Even from the result obtained in 
these two epochs when NGC 4151 is relatively bright, 



the slope variation of the Baldwin effect is still evident 
(see Fig. 4). Using the HST/STIS spectra, we also inves- 
tigate the difference in the spectral fitting with our model 
and the model consisting of three Gaussian components to 
represent three low-ionization absorption lines at around 
1526 A, 1533 A, and 1577 A, respectively. The difference is 
about 2% when the UV continuum flux is larger (F ° = 

° y 1440A 

2.4 x 10~ 13 erg s" 1 cm^A" 1 on JD2450855) and about 
13% when the UV continuum flux is smaller (F ° = 

v 1440A 

1.4 x 10~ 14 erg s" 1 cm^A" 1 on JD2451334). Although 
the equivalent widths of these low-ionization absorption 
lines become greater when the source is fainter, we see 
that the simplification of our spectral fitting does not lead 
to significant underestimation of the total C iv flux even 
in the faint state of NGC 4151. Together with the similar 
trend in the Baldwin effect found for Fairall 9, which does 
not show absorptions in Civ, and for NGC 5548, which 
docs show some narrow absorption features in C IV but 
not in Hp, we believe that the non-constant slope of the 
Civ Baldwin effect in NGC 4151 should not be driven 
mainly by the absorption effect. However, more quantita- 
tive studies of the effects of complex absorption features 
on the Baldwin effect of NGC 4151 are still needed. 

4. Conclusion and discussion 

With the long-term UV spectral data of IUE, HUT, & 
HST, we found the non-constant slope in the intrinsic 
Civ Baldwin effect of Seyfert 1 galaxy NGC 4151. The 
trend of the slope change with the UV continuum flux 
variation is similar to those found for two other Seyfert Is, 
Fairall 9 and NGC 5548, suggesting that the non-constant 
slope may be usual for AGNs. The physical origin of such 
a non-constant slope is probably related to the different 
response of the broad-line emission to the continuum vari- 
ations at different luminosity levels (Korista et al. 1998). 

From the theoretical point of view, an accretion disk 
at different luminosity levels (corresponding to different 
accretion rates) should have different accretion modes. 
Probably the AGNs accrete in a radiatively inefficient ac- 
cretion flow (see Narayan, Mahadevan & Quataert 1998 
for a review) at a lower accretion rate and in a standard 
optically thick disk (Shakura & Sunyaev 1973) or a slim 
disk (Abramowicz et al. 1988) at a higher accretion rate. 
Different accretion modes radiate differently, producing 
a relatively hard spectrum at a lower accretion rate and 
a soft spectrum at a higher accretion rate. Because the 
ionizing continuum of AGNs mainly comes from the ra- 
diation of the accretion disk, the difference in the disk- 
emitting spectrum may lead to the different response of 
the emission line to the ionizing continuum through the 
photoionization process. For NGC 4151, the lowest and 
highest continuum fluxes at 1440 A in the 1978-2002 ob- 
servational period are about 3 x 10 _15 erg s _1 cm _2 A _1 
and 6 x 10~ 13 crg s _1 cm _2 A _1 , which correspond to lu- 
minosities at 1440 A of 10 41 erg s _1 and 2 x 10 43 erg s _1 , 
respectively. If we adopt the black-hole mass for NGC 
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4151 as 1.33 x 1O 7 M (Peterson et al. 2004) and assume 
f„ oc v~ - 5 in the UV/optical band and Lh n i ~ 9L 

jv if ooi 5100 A 

for deriving the bolometric luminosity (Kaspi et al. 2000), 
we can estimate the Eddington ratio {L bo i/L Edd , usually 
taken as a measure of the dimcnsionless accretion rate) in 
the lowest and highest flux states of NGC 4151 as 0.001 
and 0.2, respectively. Clearly, the dimensionless accretion 
rate of NGC 4151 varies more than two orders of mag- 
nitude. Because the critical dimensionless accretion rate 
between a radiatively efficient accretion flow and a radia- 
tively inefficient one is about 0.01 (Narayan et al. 1998), 
most probably NGC 4151 accretes in a radiatively in- 
efficient accretion flow at the lower flux state but in a 
radiatively efficient one at the higher flux state. Such a 
change of accretion modes can not only produce different 
spectral energy distributions in the ionizing continuum, 
but can also produce different non-linear responses of the 
C iv emission line and then lead to the non-constant slope 
of the Baldwin effect in the different luminosity levels. 

Previous studies have revealed the importance of the 
light-travel time effect on the Baldwin effect (Krolik et 
al. 1991; Pogge & Peterson 1992; Peterson et al. 2002). 
Correction of such an effect can substantially reduce the 
scatters of the intrinsic Baldwin effect. However, the time 
lag for NGC 4151 between the Civ emission line and the 
UV continuum flux is not well-determined. Clavel et al. 
(1990) obtained a delay of 3.2 ± 3 days using the two- 
month data with a mean interval of 3.4 days of the IUE 
campaign in 1988-1989 when NGC 4151 was relatively 
faint. Crenshaw et al. (1996) failed to detect any correla- 
tion between the C iv line and UV continuum flux from a 
short but intensive IUE campaign in December 1993 when 
NGC 4151 was relatively bright. We also used the whole 
IUE data set in 1978-1996 and estimated the time lag 
between the total C iv line flux and the continuum flux 
at 1440 A with an interpolation cross-correlation function 
(ICCF) method developed by Gaskell & Sparke (1986), 
Gaskell & Peterson (1987), and White & Peterson (1994). 
We found that the time lag is about 1.9 days. We tried to 
use such a time lag to correct the light-travel time effect 
for NGC 4151 and found that our result does not change. 
This is mainly because the mean interval of the whole UV 
data set in 1978-2002 is about 15 days, which is signifi- 
cantly longer than the estimated time lag. 

In conclusion, we find the non-constant slope of the 
intrinsic Baldwin effect in NGC 4151 and suggest that 
such a non-constant slope may not be unusual for AGNs. 
Obviously, more studies on more AGNs are needed to con- 
firm our result. This will require intensive, long-term, and 
high-resolution UV spectral monitoring with the current 
and future space UV observatories on some strongly vari- 
able AGNs. We expect that these future studies will help 
us to extract more accurate emission-line and continuum 
fluxes and establish a more reliable line-continuum rela- 
tion for AGNs. 
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